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Abstract 

The transverse- momentum {Qt) distribution of low-mass Drell-Yan pairs 
is calculated in QCD perturbation theory with all-order resummation of 
as(as ln(Q^/Q^))" type terms. We demonstrate that the rapidity distribu- 
tion of low-mass Drell-Yan pairs at large-enough transverse momentum is an 
advantageous source of constraints on the gluon distribution and its nuclear 
dependence. We argue that low-mass Drell-Yan pairs in the forward region 
provide a good and clean probe of small-x gluons at collider energies. 
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I. INTRODUCTION 



Precise knowledge of the gluon distribution at small x is critical for reliable predictions of 
hard processes in hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions at collider 
energies. Although tremendous effort has been invested in the study of the gluon distribution 
of a free hadron, there are still large uncertainties in the distribution in both the large and 
small X regions [1-4] . Because of the limited kinematic range that can be probed by fixed 
target arrangements, we have not had enough data to shed light on the gluon distribution in a 
nucleus. Thus, information on the gluon distribution in a nucleus is far from complete [5-16]. 
With better data from coUider facihties, hke the Relativistic Heavy Ion CoUider (RHIC) and 
the future Large Hadron Collider (LHC), and also from the Tevatron at Fermilab, we expect 
much better information on the gluon distribution in free hadrons as well as in nuclei. In 
this paper, we demonstrate that the rapidity distribution of low-mass Drell-Yan pairs at 
large-enough transverse momentum is an advantageous source of constraints on the gluon 
distribution and its nuclear dependence. 

A good process from the point of view of probing the gluon distribution must satisfy 
the following criteria: (i) it has to be reliably calculable within the factorization frame work 
of perturbative Quantum Chromodynamics (pQCD); (ii) its production cross section must 
be dominated by gluon-initiated sub-processes; and (iii) it should have sufficient production 
rate for the process to be observed. For probing the gluon distribution in a nucleus, we need 
to add another criterion: (iv) the process should have relatively small initial- and final-state 
medium effects because of the twist-2 nature of the gluon distribution. 

It was first pointed out by Berger et al. [17] that the transverse-momentum distribution 
of Drell-Yan (massive) lepton pairs produced in hadronic collision is an advantageous source 
of constraints on the gluon distribution because of the dominance of gluon initiated subpro- 
cesses. These authors also emphasized that this process is free from the complications of 
photon isolation that beset studies of prompt photon production [17,18]. Due to the "Drell- 
Yan factor" in the cross section (the first factor on the right-hand side of Eq. (1) below). 
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it is useful to focus on low- mass Drell-Yan pairs to satisfy condition (iii). This, however, 
necessitates a recently-introduced resummation procedure [19]. 

In this paper we study the rapidity and transverse-momentum distribution of low-mass 
Drell-Yan pairs at RHIC, Tevatron and LHC energies, and the nuclear dependence of these 
distributions. (A short account of some aspects of the present work was given earher [20].) 
In Sec. II the theoretical framework of our calculation is briefly reviewed. In Sec. Ill, we 
show that if the transverse momentum Qt is larger than the invariant mass Q of the pair, 
the gluon initiated subprocesses give more than 80% of the low-mass Drell-Yan cross section 
at these colliders. We present our calculations for the cross sections of massive lepton-pair 
production at these coUider energies in Sec. IV, and demonstrate that low-mass Drell-Yan 
lepton pairs are measurable. We also show that the gluon dominance does not reduce for 
all relevant rapidities. In Sec. V, we study the range of parton momentum fraction x that 
can be probed at these coUiders while we vary the rapidities. We find that the forward 
region in rapidity is an excellent place to probe small-x gluons. In Sec. VI, we calculate the 
nuclear modification ratio, Rab, for dAu, AuAu, pPb and PbPb collisions. We investigate 
the effect of isospin for heavy ion beams, and the effect of nuclear modifications on parton 
distribution functions. We find that isospin plays an important role and study its impact 
on the Rab ratio. We also observe an interesting dependence and constraint on the effective 
nuclear parton distributions. We point out that the forward suppression of Drell-Yan pair 
production in nucleus-nucleus collisions at LHC energies is expected to show a very different 
behavior from the pattern at RHIC. Our conclusions are given in Sec. VII. 



II. CALCULATIONAL FRAMEWORK 

The massive lepton pair of the Drell-Yan process is produced via the decay of an inter- 
mediate virtual photon, 7*. Within the context of pQCD, the Drell-Yan cross section in a 
collision between hadrons A and B, A{Pa) + B[Pb) — > 7*(^ ll{Q)) + can be expressed 
in terms of the cross section for production of an unpolarized virtual photon of the same 
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invariant mass [17], 

daAB-^£+i-{Q)X _ ( en em \ da AB-^-i*{Q)X 



dQ'^ dQl dy \ SttQ^ j dQl dy 



(1) 



with 



da 



J dXi(f)a/A{xi, Ijl) j dX206/B(x2,//) 



dQj' dy 

x^^^^f^{x„x2,Q,QT,y;f^) , (2) 



where Tiafi runs over all parton flavors; the variables Q, Qt, and y are the invariant mass, 
transverse momentum, and rapidity of the pair, respectively; X stands for an inclusive sum 
over final states that recoil against the virtual photon. In Eq. (2), the (pa/A and (pb/B are par- 
ton distribution functions, /i represents both the renormalization and the factorization scales 
(taken to be equal), and the daah-*'y*{q)x/ dQ\dy are short-distance hard parts calculated 
perturbatively order-by-order in powers of as- In Eq. (1), an integration has been performed 
over the angular distribution in the lepton-pair rest frame. Because of the Drell-Yan factor, 
OiemI {^T^Q^) in Eq. (1), the Drell-Yan cross section suffers from a low production rate [19]. 

To increase the production rate, it is desired to measure low-mass lepton pairs. On 
the other hand, the perturbatively calculated partonic hard parts, daab-,j*{Q)x/dQ'^dy in 
Eq. (2), have q;s(q;s ln((5|,/(5^))"' type large logarithms if Q"^ <^ Q'^, and the process has a 
potential of large background from open charm (or bottom) decay. Recently, it was shown 
[19] that the perturbatively factorized Drell-Yan cross section in Eq. (2) can be systematically 
re-organized into a new factorized form so that the as{as ln{Q'^ / Q'^))'^ type large logarithms 
are resummed to all orders in a^: 

daAB^-y {Q)X j dXi(j)a/A{xi,lJi) j rfX206/B(x2, /i) 



dQj. dy 



a,b 

^\ — "!^2 f''^ (^1' ^2, Q, Qt, I/; ij) (3) 



dQjp dy 

^da^2^,x Q 
dQ% dy z 



where 'Ea,b and Ec run over all parton flavors. In Eq. (3), da'^^^^ ^^q^-^ / dQ'^ dy and 
dd'j^l^cX / dQxc dy, ai'e rc-organized perturbatively calculated short-distance hard parts with 
the superscripts (Dir) and (F) indicating the direct and fragmentation contributions, re- 
spectively; the Dc-»7*(Q)x are virtual photon fragmentation functions [21], which include 
all order resummation of the as{asln{Q'^/Q'^))'^ type large logarithms. The quantity /xr- 
represents the fragmentation scale. The four-vector is defined to be but with Q'^ set 
to be zero. The CTEQ5M parton distribution functions (PDFs) are used in our calculations 
reported here, with all scales set to \JQt + Q^- This combination of the physical scales in 
the problem appears very natural and has appropriate limits when <C and when 

<^ Qt- We are primarily interested in the latter case in the present work. 

As pointed out in Ref. [19], the resummation of the ^^(tts ln((5f^/Q^))'* type large log- 
arithms is to re-organize the perturbative short-distance hard part, daab^'y*{(:})x/dQ^dy 
in Eq. (2) into a "direct", da^^^^^g^^/dQ'^ dy, and a "fragmentation" contribution, 
do'^abLcX / dpTc dy 1 (Eq. (3)). It is important to realize that the direct and fragmentation 
contributions each have their separate perturbative expansions, recasting the single expan- 
sion of the conventional QCD factorization approach into two series [19]. A significant 
advantage of this re-organization is that the direct and fragmentation expansions are eval- 
uated at a single hard scale and are free of large logarithms, and consequently, much more 
stable perturbatively. 



III. GLUON-INITIATED CONTRIBUTION 

To study the relative size of gluon-initiated contributions, we define the ratio 

_ daAB-^r {Q)x (gluon-initiated) / daAB^riQ)x / 
'~ dQldy I dQ\dy ' ^ ^ 

The numerator includes the contributions from all partonic sub-processes with at least one 
initial-state gluon, and the denominator includes all sub-processes. In Fig. 1(a), we show 
Rg as a function of Qt at fixed y, while Fig. 1(b) displays Rg as a function of y at fixed Qt 
at Q = 2 GeV for proton-proton collision at RHIC energy. 
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FIG. 1. The ratio defined in Eq. (4) with Q = 2 GeV at ^/s = 200 GeV. Part (a): Rg 
as a function of Qt for different rapidities; part (b): Rg as a function of rapidity y for different 
transverse momenta. 

Figures 2 and 3 contain similar information for proton-antiproton collision at Tevatron 
energies, ^/s — 630 GeV and ^/s — 1960 GeV, respectively. 

The Qt dependence of Rg shows a rapid increase up to Qt ~ Q, followed by a plateau 
with a value at around (or above) 0.8 at all energies. This confirms that, for sufficiently 
large Qt, gluon-initiated sub-processes dominate the Drell-Yan cross section and thus low- 
mass Drell-Yan lepton-pair production at large transverse momentum is an excellent source 
of information on the gluon distribution at RHIC, Tevatron and LHC energies [17]. The 
fall-off of Rg towards very large Qt is related to the reduction of phase space and to the 
fact that the cross sections arc evaluated at larger values of the parton momentum fractions. 
As a function of rapidity, Rg increases slightly as one moves in the forward direction (in 
particular at RHIC energy, see parts (b) of the Figures) . This already calls attention to the 
forward region as a good kinematical domain to test gluon distributions. An example of 
the energy dependence of Rg is displayed in Fig. 4. We see that Rg increases slightly as ^/s 
increases, but the dependence is not too strong. (Note the logarithmic scale for ^/s.) 
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FIG. 2. The ratio Rg defined in Eq. (4) with Q = 2 GeV at ^/s = 630 GeV. Part (a): Rg 
as a function of Qt for different rapidities; part (b): i?^ as a function of rapidity y for different 
transverse momenta. 
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FIG. 3. The ratio Rg defined in Eq. (4) with Q = 2 GeV at = 1.96 TeV. Part (a): Rg 
as a function of Qt for different rapidities; part (b): Rg as a function of rapidity y for different 
transverse momenta. 



7 



1 

0.95 

0.9 
0.85 

0.8 
0.75 

0.7 
0.65 

0.6 



Q=2 GeV 



10 




y=0, Qt=4 GeV 
y=3, Qt=4 GeV 



10" 



VS (GeV) 



FIG. 4. The cm. energy dependence of the ratio Rg defined in Eq. (4) for Q =2 GeV, Qt = 
4 GeV, y = (solid) and y = 3 (dashed). Note the logarithmic energy scale. 



IV. THE DRELL-YAN CROSS SECTION AT RHIC AND TEVATRON ENERGIES 

Recalling criterion (iii) in the Introduction, we now turn to the production rate of low- 
mass Drell-Yan pairs at RHIC and Tevatron. We present the Drell-Yan cross section at large 
transverse momentum with all-order resummation. For Drell-Yan pair production, as seen 
from Eq. (1), there is a phase space penalty associated with the finite mass of the virtual 
photon, and the Drell-Yan factor, aem/(37rQ^) < lO'^/Q^, renders the production rate for 
massive lepton pairs small at large values of Q. In addition, the spectra drop rapidly with 
increasing Qt- In order to enhance the Drell-Yan cross section while keeping the dominance 
of the gluon initiated sub-processes, it is useful to study lepton pairs with low invariant mass 
and relatively large transverse momentum [19]. With the large transverse momentum Qt 
setting the hard scale of the collision, the invariant mass of the virtual photon Q can be 
small, as long as the process can be identified experimentally, and Q ^ Aqcd- For example, 
the cross section for Drell-Yan production was measured by the CERN UAl Collaboration 
[22] for virtual photon mass Q e [2m^, 2.5] GeV. 

Figure 5 shows the all-order resummed result for the Drell-Yan cross section as a function 
of Qt at the RHIC energy of ^/s = 200 GeV and rapidity y = for two values of the mass, 
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(a) Q — 1 GeV and (h) Q — 2 GeV. The cross section increases by about a factor 10 when the 
Drell-Yan mass decreases from 2 GeV to 1 GeV. It might still be a challenge to measure the 
low-mass Drell-Yan production with the production rate shown in Fig. 5. At the LHC, both 
the collision energy and luminosity are significantly improved. Our calculation shows that 
the production rate at the LHC is sufficiently large for being measured. (We also calculated 
Rg for (5 = 1 GeV. When Q decreases to 1 GeV, we find that Rg actually increases slightly 
in the Qt ~ 1 GeV region and it is similar to the Q — 2 GeV case in the high Qt region.) 



(V 

o 
\ 

a. 

<x 

O 
x> 

o 

b 



10 
10 
10 

1 

10 ' 
10 " 
10 " 
10 













a'"' 


r \ 
r \ 


Q=1 GeV 


: (o) 




r y=0 




: VS=200 GeV 




1 





10 15 



20 

QT(GeV) 



25 



10"^ 
10^ 
10 

1 



0) 

o 
\ 

a. 

X) 

o 

CM 

-1 

^ 10 

b 

-O -2 

10 



10 



10 



10 













a'*' 




Q=2 GeV 


: (b) 




r y=0 




: VS=200 GeV 









5 10 15 20 25 

QT(GeV) 



FIG. 5. Drell-Yan cross section as a function of Qt at RHIC energy -s/s = 200 GeV and rapidity 
y = 0, for (a) mass Q = I GeV, and (b) Q = 2 GeV. The top (solid) lines represent the total cross 
sections; leading-order (LO, dashed) and higher-order (HO, dotted) contributions are also shown. 



Figure 6 displays similar results for the Tevatron energies of -^/s = 1.96 TeV ((a) and 
(b)) and ^/s — 630 GeV (c). For the larger energy, we show results with (a) Q = 5 GeV 
and (b) Q — 2 GeV. The cross section is about 20 times larger for Q = 2 GeV than for 
Q = 5 GeV at any fixed Qt- 

In Figs. 5 and 6 we separately show the leading-order (LO, dashed), and higher-order 
(HO, dotted) contributions. At high Qt the LO and HO contributions are roughly compa- 
rable. The total cross section can also be broken up into gluon-initiated, quark-quark, and 



resummed contributions. Results on this composition are published elsewhere [20]. 
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FIG. 6. Drell-Yan cross section at rapidity y = as a function of Qt at Tevatron energies: (a) 
^/s = lM TeV and Q = 5 GeV; (b) ^/s = 1.96 TeV and Q = 2 GeV; and (c) ^/s = 630 GeV and 
Q = 2 GeV. The top (solid) lines represent the total cross sections; leading-order (LO, dashed) 
and higher-order (HO, dotted) contributions are also shown. 



Figures 7 and 8 present the rapidity dependence of the cross sections at RHIC and 
Tevatron energies at various values of Q and Qt- It can be seen that while there are order- 
of-magnitude differences in the cross section between different values of Q and Qt, the 
rapidity distributions remain flat upto high absolute values of y. This lends further support 
to the idea of exploring forward rapidities. As we will see in Section V, the forward region 
is an excellent place to probe small-x gluons. 
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FIG. 7. Drell-Yan cross section as a function of y at RHIC energy y/s = 200 GeV with 
Qt = 1 GeV (solid) and Qt = 2 GeV (dashed) for (a) Q = 1 GeV and (b) Q = 2 GeV. 
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FIG. 8. Drell-Yan cross section as a function of rapidity y for Q = 2 GeV and Qt = 3 GeV 
(solid), 7 GeV (dashed), and 9 GeV (dotted) at Tevatron energies: (a) ^/s = 630 GeV, and (b) 
^/s = 1960 GeV. 
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V. PARTON MOMENTUM FRACTION PROBED 



It is very important to estimate the region of x in the gluon distribution probed by 
the low-mass Drell-Yan process. Let us use Xi to refer to partons in one of the beams 
(say the "yellow" beam in the RHIC color code), and X2 for partons in the other ("blue") 
beam. The Xi and x^. integrations of the Dell-Yan cross section run from the appropriate 
minimum values (given by the kinematics) to 1. At central rapidity y = 0, the integration 
is symmetric in xi and X2- We are most interested in small x physics. We expect a large 
fraction of the cross section to come from rather small values of x^ in the forward region 
from the perspective of the yellow beam, where xi is not small. To quantify our statement 
with respect to X2, we introduce a cutoff Xcut to limit the X2 integration to the interval Xmin 
to Xcut for partons in the blue beam. Let us define the ratio 



Due to the symmetry in x\ and X2-, at central rapidities it is enough to examine the shape of 
to establish which region dominates the x integration — the region of x low-mass Drell- 
Yan data can provide precise information about. In Fig. 9(a) we plot Rx^ at RHIC, as a 
function of Xc^t for Q = 1 GeV for different rapidities and transverse momenta. Figure 9(a) 
shows that in the central rapidity region X2 ~ [10~^, 10~^] dominates the integration for 
both Qt = 1 GeV and Qt = 2, GeV when y = at RHIC. 

Figure 9(a) also provides important information about the forward region (from the 
perspective of Xi). At y = 3, 90% of the cross-section is given by X2 < 0.01, which is exactly 
the shadowing region of the nuclear parton distribution function. This region dominates 
the integration for both Qt — 1 GeV and Qt — 3 GeV when y = 3. Our calculation for 
the rapidity distribution of the total cross section (see Fig. 7) shows that the cross section 
does not start to drop until y ~ 3. In other words, the production rate in the large-rapidity 
region is not significantly smaller than in the central region, as long as \y\ < 3. While, as 
indicated by Fig. 9(a), the typical X2 is much smaller in the forward region than at |/ = 0, 
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the situation is different for xi. In order to evaluate the important region of xi for Drell-Yan 
production, we define the analogous ratio (compare to Eq. (5)) 




In Fig. 9(b), we plot Rx^ for forward rapidity y = 3 and Q =1 GeV at Qt = 1 GeV 
(solid) and Qt — 3 GeV (dashed). For Qt — 1 GeV, the dominant region is Xi ~ [0.1, 0.3] 
and for Qt = 3 GeV, the dominant region is Xi ~ [0.3, 0.7]. (The ratio R^^ ior y = is of 
course identical to Rx2 for y = due to the symmetry, and thus can be seen in Fig. 9(a).) As 
we will discuss later, the dominant X2 region plays an important role in the nuclear effects 
on Drell-Yan production in an AB nuclear collision. 
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FIG. 9. The ratios Rx2 (defined in Eq. (5)) and R^i (Eq. (6)) for Q = 1 GeV at RHIC energy, 
a/s = 200 GeV. In (a), both central rapidities: y = 0, Qt = 1 GeV (solid); y = 0, Qt = 3 GeV 
(dotted), and forward rapidities: y = 3, Qt = 1 GeV (dashed); y = 3, Qt = 3 GeV (dot-dashed) are 
shown. In (b) Rxi is displayed at forward rapidities: y = 3, Qt = 1 GeV (solid); ?/ = 3, Qt = 3 GeV 
(dashed). See text for more explanation. 

Figures 10 and 11 are similar to Fig. 9, but for the top Tevatron energy of =1.96 TeV 
and for the LHC energy of 5.5 TeV, respectively. Fig. 10(a) displays Rx^ as a function of 

Xcut for Q = 2 GeV for different rapidities and transverse momenta Qt- Fig. 10(b) shows 
Rxi for y = 3 at Qt — 3 GeV and Qt — 7 GeV. It can be seen from Fig. 10(a) that the 
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dominant region in X2 at central rapidity is X2 ~ [10~^, 10~^] at the top Tevatron energy; 
at forward rapidities X2 < 10~^ is reached. For xi in the forward region, the dominant 
contribution comes from xi ~ [4 * 10~^, 10~^], as seen in Fig. 10(b). In Fig. 11(a) we plot 
as a function of Xcut ior Q — 2 GeV for different rapidities and transverse momenta at 
the LHC. Figure 11(a) shows that in the central rapidity region, even when Qt — 10 GeV, 
X2 < 0.01, i.e. the shadowing region, dominates the contribution to the Drell-Yan cross 
section. In the forward rapidity region with y — 2.4, low-mass Drell-Yan production will 
probe partons with X2 as small as on the order of 10"^. In Fig. 11(b), we plot Rx^ for y = 2.4 
at Qt = 2 GeV and Qt = 10 GeV. The information provided here will be used in the next 
section (Section VI) to discuss the nuclear effects in the forward region at the LHC. 




FIG. 10. The ratios Rx2 (defined in Eq. (5)) and R^i (Eq. (6)) for Q = 2 GeV at Tevatron en- 
ergy, ^/s = 1.96 TeV. In (a), both central rapidities: y = 0, Qt = 3 GeV (solid); y = 0, Qt = 7 GeV 
(dotted), and forward rapidities: y = 3, Qt = 3 GeV (dashed); y = 3, Qt = 7 GeV (dot-dashed) are 
shown. In (b) Rxi is displayed at forward rapidities: y = 3, Qt = 3 GeV (solid); y = 3, Qt = 7 GeV 
(dashed). 
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FIG. 11. The ratios (a) R^^ and (b) R^^ for Q = 2 GeV at LHC ^/s = 5500 GeV. Different 
lines in (a) represent Rx2 at different y and Qt- y = 0, Qt = 2 GeV (solid); y = 0, Qt = 10 GeV 

(dotted); y = 2A,Qt = 2 GeV (dashed); y = 2.4, Qt = 10 GeV (dot-dashed). In (b) we have 
y = 2.4, Qt = 2 GeV (solid) and y = 2.4, Qt = 10 GeV (dashed). 

VI. NUCLEAR EFFECTS AT RHIC AND LHC 

In lack of nuclear effects on the hard collision, the impact-parameter integrated pro- 
duction cross section of low-mass Drell-Yan pairs in nucleus-nucleus {AB) collisions should 
scale, compared to the production in pp collisions, as the number of binary collisions. How- 
ever, there are several nuclear effects on the hard collision in a heavy-ion reaction. For high 
transverse-momentum Drell-Yan production, two kinds of nuclear effects are expected to 
play an important role: isospin effects and the modification of the parton distribution func- 
tion in the nucleus (shadowing effect). In this Section, we calculate a "nuclear modification 
factor" Rab at RHIC and LHC, and discuss isospin and shadowing effects in Drell-Yan pro- 
duction. We consider different rapidities, and for the RHIC program, compare AuAu and 
dAu collisions. Isospin effects were recently also emphasized in dAu scattering at forward 
rapidities in Ref. [23]. 

Similarly to the "nuclear modification factor" in hadron production, we define a ratio. 
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Rab{Qt) = 



da^'^\ZA/A,ZB/B) I da 



(7) 



dQ'^dy I dQ^dy ' 

where Za and Zb arc the atomic numbers and A and B arc the mass numbers of the 
coUiding nuclei, and the cross section da'^^'^\ZA/A, Z b / B) / dQ'^.dy uses nuclear PDFs defined 
in Eq. (8) below, while da/dQ\dy is the pp cross section defined in Eq. (2). Thus, Rab 
incorporates the effects of both shadowing and nuclear composition (isospin effect). 
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FIG. 12. The ratio Rab defined in Eq. (7) for (a) Q = 1 GeV and (b) Q = 2 GeV. Different 
fines in the figure represent Rab at different values of y for different collisions. 



In Fig. 12, we plot Rab as a function of Qt at RHIC in dAu and AuAu collisions at both 
central rapidity {y = 0) and forward rapidity {y = 3) for {a) Q = 1 GeV and {h) Q = 2 GeV. 
In this work we use the EKS parameterization to incorporate the shadowing effect in the 
nuclear PDF [8]. Numerical differences from using other nuclear PDF parameterizations 
[5,10-13,15] could be estimated from the ranges of momentum fraction x probed at different 
collision energies. Fig. 12 shows that in the central rapidity region, RdAu is larger than RauAu 
in the small Qt region, and RauAu becomes larger than one and larger than RdAu when Qt 
is larger than ki 3 GeV. This seems to be surprising at first glance. Fig. 12 also shows strong 
suppression of Drell-Yan pair production in both dAu and AuAu collisions in the forward 
region. Furthermore, RauAu increases steeply at large Qt > 6 GeV when y = 3. In order to 
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disentangle shadowing and isospin effects, botfi incorporated in Rab as defined in Eq.(7), 
we separate tfiese dependences and investigate isospin and sfiadowing effects separately in 
the following discussion. 

To take into account the neutron-proton composition of nuclei one normally uses 



for the nuclear PDF, where ^^/a proton PDF and ^^/a PDF of the neutron. 

By artificially setting Za = A and Zb = B m our calculation, we consider nuclei made of 
protons only. Calculating the analogue of (7) for these hypothetical nuclei we define -R^'^, 
which is free of isospin effects and contains shadowing only (pure shadowing). The ratio 
Rab shown at |/ = in Fig. 13 by the dotted lines for (a) dAu and (b) AuAu collisions. 

It can be seen from Fig. 13 that in the central rapidity region, when Qt > 3 GeV, Ra^ 
is larger than 1 in both AuAu and dAu collisions. This means that anti-shadowing in the 
PDF start to be important at Qt ~ 3 GeV. The result is consistent with what was learned 
about the dominant region of x earlier. Recall from Fig. 9(a) that the dominant region of x 
for partons from both beams is about [10~^, 10~^] for Qt — 3 GeV and y — Q. An inspection 
of the EKS nuclear PDFs [8] informs that x ~ 10"^ is the region where anti-shadowing 
starts to take over shadowing. In AuAu collisions, partons in both beams experience the 
same shadowing (anti-shadowing) effects at central rapidity. In dAu collisions, only the 
partons in the Au nucleus have significant shadowing (anti-shadowing) effects. Therefore 
pure shadowing (anti-shadowing) effects are stronger in the AuAu collisions than in the dAu 
collisions as seen in Fig. 13. 

Now let us also investigate the pure isospin effects on low-mass Drell-Yan production. 
Isospin effects originate in the difference between the parton distributions of neutrons and 
protons. Due to this difference, the production cross section of Drell-Yan pairs in proton- 
neutron and neutron-neutron interactions differs from that in pp collisions. To study isospin 
effects, let us introduce 



A 



a/A 



n 



a/A 



(8) 
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where the numerator now does not include shadowing, i.e. standard CTEQ5M PDFs are 
used without shadowing. Examining R^b completes our program of separating the isospin 
and shadowing dependences of Eq. (7). 
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FIG. 13. The pure shadowing (dotted) and isospin (dashed) components of ratio Rab at y = 
and Q =1 GeV at RHIC for (a) dAu and (b) AuAu collisions. See text for more explanation. 

Figure 13 also displays (dashed hne) as a function of Qt in (a) dAu and (b) AuAu 
collisions at y = 0. We see that R^X°b < 1 in the whole transverse momentum region 
[1, 10] GeV; isospin effects suppress the production of high transverse- momentum low-mass 
Drell-Yan pairs at RHIC. The suppression is getting slightly stronger with increasing Qt- 
Since the Compton type (qg) processes are the dominant sub-processes in low-mass Drell-Yan 
production at high transverse momentum, we next consider the isospin effects in these sub- 
processes to help us understand the suppression result shown in Fig. 13. The contribution 
of different flavor quarks is proportional to the square of their fractional charges (charge 
factors): 4/9 for u quarks and 1/9 for d quarks. Let us first look at proton-proton [pp) 
and proton-neutron [pn) collisions. In the calculation of the pp Drell-Yan production cross 
section, the PDF part for the Compton-typc processes, after including the charge factors, is 

proportional to $g(4$„/9-|-$d/9-|-4$^/9-|-$j/9-| ), plus an "exchange" term, where the 

gluon comes from the other proton. For pn, when the gluon comes from the proton, we have 
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$p($u/9 + 4$d/9 + $s/9 + 4$j/9H ) (where • • • represents the contributions which are the 

same for pp and pn), plus an exchange term. In Fig. 14, we plot $^(0;, /i)/$„(a;, /x) (soUd), 
^^{x, iJ,)/^u{x, 11) (dashed), and the ratio of the above two charge-weighted expressions 
(dotted) as functions of x at the scale //^ = 10 GeV^. For not-too-small x, in particular 
in the region of [10""^, 10~^], which is the dominant region of parton x, valence quarks are 
more important than sea quarks. Since $^(0;, /x)/$u(a;, /x) (solid) is itself smaller than one, 
the production cross section in pn is smaller than in pp collisions. The dotted line in Fig. 14 
illustrates the effect. Using a similar argument, one finds that the isospin effect in nn 
collisions is a little stronger than in pn collisions. Returning to nuclear reactions, where a 
large fraction of the nucleon-nucleon collisions is not pp, but pn and nn, we conclude that 
due to the charge factor and the fact that < in the x region of interest, i?^^ is smaller 
than 1. When the transverse momentum increases, x in the PDFs increases. Since the ratio 
of becomes smaller at larger x, R^Xb decreases as Qt increases. Since the change of 

the proton/neutron ratio from d to Au is small for these purposes, B^Xb ^^^7 similar in 
dAu to that of AuAu, as Fig. 13 demonstrates. 
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To summarize the above discussion, at central rapidity y — 0, when Qt > 3 GeV, low- 
mass Drell-Yan pair production at RHIC suffers anti-shadowing. The anti-shadowing effects 
are stronger in AuAu collisions than in dAu collisions; isospin effects suppress the production 
of low-mass Drell-Yan pairs and the suppression increase with Qt] isospin effects are similar 
in dAu and AuAu coUisions. The combined effects of isospin and shadowing lead to a larger 
Rab in AuAu than in dAu at RHIC energy. 

Figure 15 displays results for Rab (solid), R^Xb (dashed), and R^J^^ (dotted) in the forward 
region for (a) dAu and (b) AuAu coUisions. It is clear that the steep rise of Rab at large 
transverse momentum in AuAu collisions is a consequence of the effects of the modification 
of the nuclear PDF in the "yellow" beam (for which y — 3 designates the forward region); 
correspondingly xi enters the Fermi-motion region of the nuclear PDF. The large suppression 
due to isospin effects, which reaches ~ 40% at high transverse momentum is also from the 
"yellow" beam side, where the larger xi gives a stronger u, d asymmetry. To demonstrate 
this point, we also plot RpAu in Fig- 15(a) (dot-dashed). The result RpAu ~ 1 confirms that 
the large isospin effect in dAu is from the d side. It is also important to point out the big 
difference between RpAu and RdAu, which illustrates that, in regard to certain observables, 
dAu collisions can be very different from pA collisions. 
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FIG. 15. The ratio Rab at RHIC in the forward region (y = 3) for (a) dAu and (b) AuAu 
collisions at Q =1 GeV as a function of Qt- 
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Figure 16 shows Rab as a function of Qt at the LHC in pPb and PbPb coUisions at 
both, central rapidity and forward rapidity. In pPb coUisions the Drell-Yan pair production 
is suppressed more in the forward region, as expected. However, it is surprising that in 
PbPb coUisions, the DreU-Yan pair production in the forward region is actuaUy sUghtly less 
suppressed than in the central rapidity region, up to Qt ~ 13 GeV. This behavior is very 
different from what was found at RHIC. It can be explained by the very different values 
of xi and X2- Prom the central rapidity region to the forward rapidity region, xi increases 
and X2 decreases. Prom Pig. 11, approximately, X2 changes from on the order of 10"^ to on 
the order of 10"'' if one wants to produce a low-mass Drell-Yan pair not at y = but at 
y — 2.4. At the same time, Xi changes from on the order of 10~^ to on the order of 10~^. 
In the EKS parameterization, the shadowing factors are almost constant in [10~^, 10~^]. 
However, shadowing effects decrease much in the region [10~^, 10~^]. This leads to the slight 
increase of Drell-Yan pair production in the forward region compared to the central rapidity 
region. Gluon saturation (color glass condensate) effects [20,25,26] may make the increase 
even larger. 
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FIG. 16. The ratios Rab at the LHC at Q =2 GeV for pPb at y = (dotted), pPb at y = 2.4 
(dash-dot), PbPb at y = (soUd), and PbPb at y = 2.4 (dashed) as a function of Qt- 




21 



VII. SUMMARY AND CONCLUSION 



In summary, low-mass Drell-Yan dilepton production is a potentially clean probe of small- 
X gluons, without strong final-state interactions. When Qt > 2 GeV, the gluon initiated 
sub-processes contribute more than 80% to the cross section. Unfortunately, low-mass Drell- 
Yan dilepton production suffers from low production rate at RHIC due to the Drell-Yan 
factor. At the Tevatron and LHC energies, high-Qr low-mass Drell-Yan production is an 
excellent probe of the gluon distributions. 

We have shown that the production rate for low-mass Drell-Yan pairs at Tevatron and 
LHC energies does not reduce for a wide range of rapidity, and the forward region is very 
sensitive to small-x gluons. Therefore, the rapidity distribution of low-mass Drell-Yan pairs 
at large enough transverse momentum could be a good and clean probe of small-x gluons at 
collider energies. 

Because of the lack of final-state interactions and the relatively weak initial-state nuclear- 
dependent power corrections at high Qt, the rapidity and transverse momentum distribution 
of low- mass Drell-Yan pairs in hadron- nucleus and nucleus-nucleus collisions should be a very 
good probe of the gluon distribution in a nucleus. By comparing the nuclear modification 
factor of pAu and dAu, we demonstrated quantitatively that the isospin effect is very im- 
portant in extracting nuclear parton distributions. We have shown that low-mass Drell-Yan 
pairs in the forward region at Tevatron and LHC energies provide very sensitive information 
on the gluon distribution and its nuclear dependence for x less than 10~^. With such sensi- 
tivity on gluons, low mass Drell-Yan dilepton can be a good probe of the novel phenomena of 
gluon saturation or color glass condensate in a large nucleus. We believe that in view of the 
importance of small-x gluons, the ability of low-mass Drell-Yan pairs to provide information 
on these gluon distributions as demonstrated in the present paper warrants the mounting 
of experimental efforts to measure low-mass Drell-Yan pairs at colliders. 
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